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ABSTRACT 

Aims. We measured radial metallicity profiles for a sample of « 50 hot, intermediate redshift galaxy clusters, selected from the XMM- 
Newton archive. 

Methods. As in our previous paper, we used background modeling rather than background subtraction, and the Cash statistic rather 
than the^ 2 . This method requires a careful characterization of all background components. We also performed montecarlo simulations 
to assess systematic effects. 

Results. The mean metallicity profile shows a peak in the center, and gently declines out to 0.2 7? [g0 - Beyond 0.2 _/? 180 the metallicity is 
^ 0.2 solar and, at variance with recently published expectations based on simulations, consistent with being flat. We find no evidence 
of profile evolution from z = 0.1 to z = 0.3. When comparing our mean profile to those obtained by recent works with BeppoSAX and 
Chandra, we find remarkable agreement over the entire radial range. 
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1. Introduction 

Clusters of galaxies are the most massive gravitationally bound 
systems in the universe. They are permeated by the hot, X-ray 
emitting intra-cluster medium (ICM), which represents the dom- 
inant baryonic component. The mean meta llicity of the ICM i s 
found to be roughly half of the solar value (Loewenstein 2004), 
therefore a substantial fraction of plasma is not of primordial 
origin. As heavy elements are only produced in stars, the pro- 
cessed material must have been ejected by cluster galaxies into 
the ICM. In rich clusters the plasma reaches temperatures of sev- 
eral 10 7 K and emits X-rays mainly via thermal bremsstrahlung. 
At such temperatures, most elements are either fully ionized or 
, in a high ionization state. The strong emission lines due to the 
transitions to the n=l level of the H-like and He-like ions of Iron 
around 6.7 keV (rest frame energy) are the most prominent fea- 
tures in X-ray spectra of hot clusters. In fact, for hot clusters as 
those belonging to our sample, the measurement of the metallic- 
ity is a measurement of the iron abundance. 

Spatially resolved analysis of metal abundance in clus- 
ters has become possible only recently, first with ASCA 
and BeppoSAX and then with Chandra and XMM- N ewton . 
These measurements (e.g. | De Gr andi & Molendi 2001; 
iDe Grandi etafl 12004 iTamura et alJ T2004; B aldi et all (20071) 
have shown that abundance gradients are common features in 
clusters of galaxies: most clusters show a peak of metallicity in 
the center and a gentle decline outward. 

In this paper we present radial metallicity profiles from the 
same sample of hot intermediate redshift clu sters analyzed in our 
previous paper dLeccardi & Molendi [2008). hereafter Paper I. 
We used the same novel data analysis technique, which is ex- 
tensively described in Paper I and summarized in Sect. [3] In 
Sect. 14.11 we show the mean metallicity profile, in Sect. 14.21 
search for an indication of the evolution of the metallicity with 



redshift, and in Sects. [431 and |4.4| compare our mean profile with 
those obtained from hydrodynamic simulations and from previ- 
ous observational works. In the Appendix lAl we describe a sys- 
tematic effect, which could affect the measurement of the metal- 
licity especially for the case of low metallicity, statistically poor 
spectra. 

Quoted confidence intervals are 68% for one interesting pa- 
rameter (i.e. AC =1), unless otherwise stated. All results are 
given assuming a ACDM cosmology with Q m = 0.3, Qa = 0.7, 
and H = 70 km s -1 Mpc~'. 



2. The sample 

We selected from the XMM-Newton^ archive a sample of hot 
(kr > 3.3 keV), intermediate redshift (0.1 < z < 0.3), and high 
galactic latitude (\b\ > 20°) clusters of galaxies. We retrieved 
from the public archive all observations of clusters satisfying the 
above selection criteria, performed before March 2005 (when the 
CCD6 of EPIC-MOS1 was switched offl) and available at the 
end of May 2007. We excluded from the sample observations 
that are highly affected by soft proton flares and observations 
of clusters that show evidence of recent and strong interactions. 
In Table 3 of Paper I we list the 48 observations that survived 
our selection criteria and report cluster physical properties. Each 
observation was performed by using THIN 1 or MEDIUM filters, 
and its total (i.e. MOS1+MOS2) exposure time is greater than 
16 ks. 



1 An ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and NASA 

2 http://xmm.esac.esa.int/external/xmm_news/items/MOS 1-CCD6/ 
index, shtml 
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3. Data analysis 

In this section we only recall the major steps of our data analy- 
sis procedure and refer to Paper I for a detailed description. We 
stress that in our procedure we used only EPIC-MOS data, be- 
cause a robust characterization of EPIC-pn background was not 
possible, mainly due to the small regions of the detector outside 
the field of view, and to the non-negligible fraction of out-of- 
time events (see Appendix B in Paper I for further details). 

3.1. Spectra preparation 

Observation data files (ODF) were retrieved from the XMM- 
Newton archive and processed in a standard way with the 
Science Analysis System (SAS) v6.1. 

The soft proton cleaning was performed using a double fil- 
tering process, first in a hard (10-12 keV) and then in a soft 
(2-5 keV) energy range. We filtered the event file according to 
PATTERN and FLAG criteria, and excluded the "bright" CCDs, i.e. 
CCD-4 and CCD-5 for MOS 1 and CCD-2 and CCD-5 for MOS2 
(see Appendix A in Paper I for a discussion). Bright point-like 
sources were detected, using a procedure based on the SAS task 
edetect_chain, and excluded from the event file. 

The Rsb indicator, i.e. the ratio between surface-brightness 
calculated inside and outside the field of view (see Eq. 1 in 
Paper I), allowed us to quantify the amount of the quiescent 
soft proton (QSP) component survived the double filtering pro- 
cess. Values of 7?sb roughly span the range from 1.0 (negligi- 
ble contamination) to 1.5 (high contamination). Since the obser- 
vation 0084230401 of Abell 267 is extremely polluted by QSP 
(^sb = 1 -8), we excluded it from the sample. 

The cluster emission is divided in 10 concentric rings 
(namely 0'-0.5', 0.5'-l', l'-1.5', 1.5'-2', 2'-2.75', 2.75'-3.5', 
3.5'-4.5', 4.5'-6', 6'-8', and 10'-12') at fixed angular radii, to 
maintain under control systematics related to the detectors (see 
the Appendices of Paper I). The most external ring was used to 
estimate background parameters only. The width of most of the 
rings is limited by the PSF of the XMM-Newton telescopes. The 
center of the rings was determined by surface-brightness isocon- 
tours at large radii and is usually, but not necessarily, coincident 
with the X-ray emission peak. We prefer that azimuthal symme- 
try be preserved at large radii, where we are interested in char- 
acterizing profiles, at the expense of central regions. For each 
instrument (i.e MOS1 and MOS2) and each ring, we accumu- 
lated a spectrum and generated an effective area (ARE). For each 
observation we generated one redistribution function (RMF) for 
MOS1 and one for MOS2. We performed a minimal grouping to 
avoid channels with no counts, as required by the Cash statistic. 

3.2. Spectral analysis 

The spectral analysis followed a substantially nevvQ approach: 
we used the background modeling, rather than the subtraction, 
and the Cash statistic rather than the x 2 - We fit spectra with 
an absorbed thermal (WABS*MEKAL in XSPEC v 11.30) plus 
background model in the 0.7-10.0 keV energy band, which rep- 
resents a good trade off between statistical quality and level of 
systematics. The details of the background model are reported in 
the Appendices of Paper I. To model the background, a careful 
characterization of all its components is mandatory. Ideally, one 








3 A somewhat similar procedure was already used by Stanfor cfet alj 
d200lh 

4 http://heasarc.nasa.gov/docs/xanadu/xspec/xspecl 1/index.html 



Fig. 1. Radial metallicity profi les for all clusters in our sa mple. 
Abundances are expressed in lAnders & Grevessel ([1989) solar 
values and radii in units of /? 180 - 



would like to estimate background parameters in the same region 
and at the same time as the source. Since this was not possible, 
we estimated background parameters in the external 10'-12' ring 
and rescaled them in the inner rings, by making reasonable as- 
sumptions on their spatial distribution tested by analyzing blank- 
field observations. 

We fit spectra in internal rings leaving the temperature and 
the normalization free to vary; the metallicity was constrained 
between +5 Z Q (see the discussion in Appendix A). The red- 
shift was constrained between ±7% of the optical measurement 
in the two innermost rings and, in the other rings, was fixed to 
the average value of the first two rings, by considering indepen- 
dently MOS 1 and MOS2 spectra. The main reason for our choice 
is to allow for EPIC calibration uncertainties, and for possible 
discrepancies between X-ray and optical derived redshift values. 
Typical shift values are on the order of 2%. Finally, we produced 
metallicity profiles for each cluster, by computing the average 
over the two MOS values. 



4. Metallicity profiles 

4.1. The mean profile 

In Fig. Q] we show radial metallicity profiles for all clusters 
in our sample. All our metallicit y measurements are relativ e 
to the solar values published by lAnders & Grevess e (1989). 
Even if these have been superseded by more r ecent val- 
ues dGrevesse & Sauvall 119981: lAsplund et ail 120051) . they al- 
low straightforward compar i son with mos t of the literature 
(iDe Grandi & Molendil 120011: iBalestra et all l2007t iBaldi et al l 
120071) . iGrevesse & Sauvall (Il998l) and lAsplund et alj (120051) 
introduced a 0.676 and 0.60 times lower iron solar abun- 
dance respectively, while other elements are substantially un- 
changed. A si mple scaling by 0.676 and 0.60 converts mea- 
sures from the Anders & Grevessel d 19891) iron abundance to 
the IGrevesse & Sauvall (1 19981) and lAsplund et al.l (120051) abun- 
dances. 

The radius is rescaled by R\sq, i.e. the radius encompassing 
a spherical density contrast of 180 with respect to the critical 
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Fig. 2. Mean metallicity pro file averaged over all cl usters. 
Abundances are expressed in lAnders & Grevessel JT989) solar 
values and radii in units of R\so- The dotted lines show the one- 
sigma scatter of the values around the average. 

Table 1. Mean metallicity profile. 



Ring" 



Metallicity" 



0.00-0.04 
0.04-0.08 
0.08-0.12 
0.12-0.18 
0.18-0.24 
0.24-0.32 
0.32-0.45 



0.450+0.010 
0.352±0.009 
0.294±0.010 
0.262±0.010 
0.230±0.013 
0.245±0.017 
0.218±0.023 



Notes: " in units of R\m\ b in solar units dAnders & Greve sse 1989). 



is concerned, we note that the size of central bins (i.e. 30") is 
comparable to the XMM-Newton PSF. As far as the recasting is 
concerned, the one we show in Fig.[2]is the result of various tri- 
als specifically aimed at reaching the best compromise between 
resolution and statistical quality; moreover, the size of central 
bins (i.e. 0.04 R\%a) corresponds to « 30" for clusters at z ~ 0.2 
wit h a temperatur e of 7 ke V. 

iBalestra et all d2007l) and iMaughan etaT] d2008l) have 
adopted a procedure that is alternative to ours to estimate mean 
metal abundances for a sample of clusters. These authors per- 
formed a simultaneous spectral fit, leaving temperature and nor- 
malization free to vary for each object and using a unique metal- 
licity value for all clusters in each redshift bin. We note that 
when modeling the background, as we do, a joint fit is infeasible, 
because of the large number of model parameters. Moreover, we 
are not aware of any detailed work that investigates the impact of 
systematic errors, possibly affecting individual measurements, 
on the final result of a joint fit. 

The mean metallicity is 0.45 Z Q in the center and decreases 
out to « 0.2 R\so\ beyond this radius the profile is consistent 
with being fiat, a fit with a constant for R > 0.2 R\$q gives 
Z = 0.23 + 0.01 Z . The profiles show a large scatter, which is 
mostly of statistical origin. In the central regions (R < 0.2 /?igo) 
we find an intrinsic scatter of 22 + 2% related to the presence 
of cool core clusters; in the outer regions (R > 0.2 R\&o) the 
intrinsic scatter is only 14 + 8% (i.e. « 0.03 Z ), the same 
order of magn itude as our systematics (s ee Appendix [At. Past 
works (e.g. |Pe Grandi & Molendil |2001|) have shown that the 
abundance profiles of cool core and non cool core clusters (see 
Sect. 4 in Paper I) differ in the central regions. We found quali- 
tatively similar results, but we choose to extensively address this 
important issue in a forthcoming paper (Leccardi et al., in prep.), 
where we also compare our results with those obtained from a 
local sample. 



density. We compute R\gQ from th e mean temperature, kTM, 
the redshift, z, as in Paper I and in lArnaud et alj (120051) : 



tfiso = 1780 



5keV 



1/2 



h(z) 1 kpc, 



and 



(1) 



where h(z) = (O m (1 + z) 3 + £2 A ) 1/2 . /? L go is a good approx- 
imation to the virial radius in an Einstein-De Sitter universe 
and has been largely used to rescale cluster radial proper- 
ties (e . g.lDe Grandi et aD2004IVikhlinin et alj|2005t iPratt et alj 
120071; iLeccardi & Molendil 120081) . The profiles are limited to 
~ 0.4 /f i8o because beyond this radius the source-to-background 
count rate ratio is too small and the measurements are unreliable 
(see Sect 3.2.2 in Paper I). 

In Table [TJ and in Fig. [2] we report the mean metallicity pro- 
file binned in units of Rmu these values have been computed by 
performing, for each new bin, a weighted average of metallicity 
values in the original bins which have a non-zero intersection 
with the new bin. The weight is the product of two components: 
one is the inverse squared errors, the other depends on the inter- 
section between the original bin and the new one. If the origi- 
nal bin is totally included into the new one, the weight is equal 
to one. If the original bin has only a partial intersection with 
the new one, the weight is the fraction of the original bin that 
belongs to the new one. Possible blurring effects associated ei- 
ther to the original binning in angular units or to the recasting in 
units of R i8o should be minimal. As far as the original binning 



4.2. Redshift evolution 

We divided clusters in our sample into two groups to investigate 
a possible profile evolution with redshift: near (distant) clusters 
are characterized by a redshift lower (greater) than 0.2. In Fig. [3] 
we show the mean metallicity profiles for the two groups. The 
single points are consistent within one- or two-sigma, except for 
the core region where the spatial resolution plays an important 
role. In our procedure we fixed the size of the central ring to 
30", which corresponds to « 0.03 Rwo for nearest {z ~ 0.1) 
and to « 0.07 Riw for most distant (z ~ 0.3) objects; thus, the 
metallicity peak for distant clusters is blurred over a larger re- 
gion. The three-sigma discrepancy in the region within 0.2 R^o 
(see shaded regions in Fig. [3]) is most likely due to a different 
fraction of cool core and non cool core clusters within the two 
subsamples (Leccardi et al., in prep.). 

In the outer regions, if the metallicity is allowed to vary be- 
tween and 5 Z as commonly done, we measure a discrep- 
ancy between two profiles of 0.061 + 0.018 Z with a more than 
three-sigma significance (see Appendix lAl. Conversely, if the 
metallicity is allowed to vary between +5 Z as in our data anal- 
ysis procedure (see Sect. 13.21 1. near and distant clusters have a 
mean metallicity of 0.248 + 0.011 Z Q and 0.219 + 0.015 Z , re- 
spectively. For this case, the discrepancy (i.e. « 1.5 sigma, see 
shaded regions in Fig. [3]) is consistent with a purely statistical 
fluctuation, and is the same order of magnitude as our system- 
atics. Summarizing, we find no evidence of metallicity profile 
evolution from z = 0.1 to z = 0.3. We stress that the system- 
atic effect described in the Appendix|A] which affects in particu- 
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Fig. 3. Mean metallicity profiles for near (z < 0.2, filled circles) 
and distant (z > 0.2, empty circles) clusters. The dark and the 
light shaded regions indicate the average metallicity within one- 
sigma uncertainties below and beyond 0.20 R^o for near and 
distant clusters, respectively. The radii have been slightly offset 
in the plot for clarity. 



lar low metallicity, statistically poor spectra, if unaccounted for, 

could cause a false detect ion of the meta l licity e volution. 

Tw o recent works (Bal estra et alj l2007t iMaughan et alj 
2008) have investigated the evolution in the iron content of the 
ICM, by analyzing data from the Chandra archive. As men- 
tioned in Sect. 14. II these authors performed a simultaneous spec- 
tral fit, leaving temperature and normalization free to vary for 
each object and using a unique metallicity value for all clus- 
ters in each redshift bin. Unfortunat ely, we cannot com pare 
our results with those obtained by Bales tra et al.l d2007l) . be- 
cause the y considered differe nt regions from cluster to cluster. 
Instead, Maug han et al.l d2008l) analyzed the region within 
(« 0.6 ^iso) with and without the core region (i.e. 0.15 /?5oo) 
and obtained respectively Z « 0.4 Z G and Z « 0.35 Z between 
z — 0.1 and z = 0.3. When analyzing roughly the same regions 
for our clusters, we obtain Z = 0.32 Z and Z = 0.26 Z respec- 
tively, with negligible uncertainties; i.e., a significantly lower 
(w 20%) mean metallicity. A possible explanation for this dis- 
crepancy may be related to different weights on the averaging 
procedure: we averaged weighting over the inverse squared er- 
rors, while iMaughan et alj d2008l) measured metal abundances 
from individual spectra extracted from the entire region of inter- 
est. 



4.3. Comparison with hydrodynamic simulations 

We compare our (hereafter LM08) mean metallicity profile with 
the one (hereafter F08) de rived from hydrod ynamic simulations 
of four relaxed clusters by Fa bian et all (2008). 

The simulations are perfor med using the hydrodynamical 
TREE-SPH code GADGET-2 dSpringelll2005l) with the imple- 
mentation of chemical enrichment by Tornato re et al.l (|2007). 
The authors used the emission-weighted definition of metallic- 
ity, with emissivity of each gas particle computed in the 0.5- 
10.0 keV energy band. In principle, for a comparison with ob- 
servational data, one should extract synthetic spectra from the 
simulated clusters and then measure the metallicity by fitting 
these spectra with a single-temperature and single-metallicity 



Fig. 4. Comparison between our observed mean profile (cir- 
cles) and the one deriv e d from hydrodynamic simulations (solid 
line) by [ Fabian et al. (2008). Abundances are expressed in 
lAnders & Grevesse ( 1989) solar values and radii in units of R\go. 



plasma model (e.g . the MEKAL mo del in XSPEC). A recent 
work presented by iRasia et al.l d2008l) showed that, at least for 
Iron, the emission-weighted estimator of the metallicity gives 
results quite close (within about 10%) to those obtained from 
the spectral-fitting analysis. 

In Fig. [4] we compare LM08 and F08 profiles. We note dif- 
ferences at both small and large radii. In the central region, the 
F08 profile is much more peaked; how ever, this difference likely 
results from two factors: namely that Fa bian et"aT] ((2008) ana- 
lyzed only relaxed clusters and that our data are limited by the 
XMM-Newton PSF At large radii, while the F08 profile shows 
a constant decrease, ours is consistent with being flat beyond 
» 0.2 i?i8o- This discrepancy could be due to issues related to 
observations, simulations, or both. More precisely, observations 
could be affected by unknown systematic effects, and simula- 
tions could underestimate the metallicity at large radii by under- 
estimating possible convective motions in the ICM, which might 
be responsible of mixing metals. 

4.4. Comparison with previous works 



We co mpare o ur results w i th thos e obtained bvlDe Grandi et ail 
(2004|) and iBaldi et al.l d2007l) . iDe Grandi et al.l (hereafter 
DM04) have analyzed a sample of 21 hot (kT > 3.5 k eV), nearby 
(z < 0.1) galaxy clusters observed with BeppoSAX. Bald i et all 
(hereafter BA07) have analyzed 12 very hot (kT > 6 keV), inter- 
mediate redshift (0.1 < z < 0.3) clusters observed with Chandra. 

Comparing results obtained from different works is not triv- 
ial; indeed, cluster physical properties, instrumental characteris- 
tics, and data analysis procedures may differ. Moreover, each 
author uses his own recipe to derive a scale radius. We have 
rescaled DM04 and BA07 profiles by using the standard cos- 
mology (see Sect. Q} and deriving the scale radius, Ri&o, as ex- 
plained in Sect. 14. It the aim is to reduce all inhomogeneities as 
much as possible. We als o c onverted DM04 abundances from 
lGrevesse& Sauvall d!998l) to lAnders & Grevessd dl989l) solar 
values. 

In Fig. [5] we compare LM08, DM04, and BA07 mean metal- 
licity profiles. Even if cluster samples cover a different redshift 
range and the instruments (i.e. XMM-Newton, BeppoSAX, and 
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Fig. 5. Mean metal licity profiles obtaine d from this work 
(LM0 8, circles) , by |Pe Grandi et alj d2004h (DM04, triangle), 
and iBaldietail d200lT (BA07 , squar es) Abundances are ex- 
pressed in lAnders & Grevessel (fl989) solar values and radii in 
units of /?i8o- The radii have been slightly offset in the plot for 
clarity. 



Chandra) present different characteristics, the mean profiles are 
remarkably similar over the entire radial range. 



5. Summary and conclusions 

We analyzed a sample of « 50 hot, intermediate redshift galaxy 
clusters to measure their radial properties, in this paper we fo- 
cused on metallicity profiles. Our main results are summarized 
as follows. 

- The mean metallicity is 0.45 Z in the center and decreases 
out to a 0.2 /?i8o; beyond 0.2 R^q the metallicity is consis- 
tent with being flat at 0.23 + 0.01 Z . 

- The profiles show a large scatter, which is mostly of statisti- 
cal origin. In the central regions the scatter (i.e. 22 + 2%) is 
also related to the presence of cool core clusters, in the outer 
regions it (i.e. 14 + 8%) is comparable to systematics. 

- There is no evidence of profile evolution from z = 0.1 to 
z = 0.3. 

- We find a mean metall icity as 20% lower than found 
by iMaughan et al.l (2008) using Chandra, when analyzing 
roughly the same regions for clusters in the same redshift 
range. 

- We point out the existence of a systematic effect, affecting in 
particular distant clusters, which, if unaccounted for, could 
cause a false detection of evolution. 

- When comparing our mean pr ofile to the one der ived from 
hydrodynamic simulations by iFabian et alj d2008l) . we find 
differences at small and large radii. In particular, while the 
profile obtained by Fabia n et ail ((2008) shows a constant de- 
crease, ours is consistent with being flat beyond » 0.2 /?igo. 

- When comparing our mean p rofile to those obtaine d by 
recent w orks with BeppoS AX dDe Grandi et al. 2004) and 
Chandra (Bal di et ai]l2007l) . we find remarkable agreement 
over the entire radial range. 

Our results have been obtained using the same novel data 
analysis technique as described in Paper I; i.e., the background 



modeling rather than the background subtraction, and the Cash 
statistic rather than the^f 2 . 
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Appendix A: Discussion of systematics 

We made use of montecarlo simulations to test the reliability 
of our metallicity measurements. The simulation procedure is 
similar to that described fo r the source-only case in Sect. 2 of 
iLeccardi & Molendil d2007l) . We considered a thermal (MEKAL 
in XSPEC) spectrum only, without a background. Input param- 
eters are 6 keV temperature, 0.25 Z metallicity, 0.2 redshift, 
and 9 x 10~ 4 (in XSPEC units) normalization. We used the 
Abell 1689 EPIC-MOS1 observation as a guideline, for produc- 
ing RMF and ARF, and for choosing typical input model pa- 
rameters. We considered the 3.5'-4.5' ring only. For each chan- 
nel, we perturbed the number of counts with a Poisson distribu- 
tion centered on the expected value, and repeated this procedure 
1000 times to obtain 1000 spectra, which simulate 1000 inde- 
pendent measurements of the source. We fit simulated spectra 
with a MEKAL model in the 0.7-10.0 keV energy band using the 
Cash statistic. Temperature, metallicity, redshift, and normaliza- 
tion were allowed to vary within their XSPEC standard ranges. 
Here we focused on the metallicity, for which the standard al- 
lowed range is between and 10 3 in solar units. For each mea- 
surement, we determined the best fit value and the one-sigma 
uncertainties. 

In Fig. IA.ll we report the frequency distribution for metallic- 
ity best fit values for a 100 ks exposure time (i.e. for spectra with 
high statistical quality); as expected, the curve is very similar to 
a narrow Gaussian peaked around the input value (i.e. 0.25 Z Q ). 
The mean, the median, and the weighted (over one-sigma uncer- 
tainties) average are all close to 0.25 Z , namely the weighted 
average is 0.246+0.002 Z . In the left panel of Fig. lAT2l we report 
the same histogram for a 20 ks exposure time (i.e. for spectra 
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Fig. A.l. Frequency distribution for metallicity (expressed in so- 
lar units) best fit values for a 100 ks exposure time. The dashed 
line indicates the input value, i.e. 0.25 Z . 
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Fig. A.2. Frequency distribution for metallicity (expressed in so- 
lar units) best fit values for a 20 ks exposure time. The dashed 
lines indicate the input value, i.e. 0.25 Z . In the left panel, the 
metallicity is forced to be positive; in the right one, it can assume 
negative values too. The two histograms only differ around zero. 



roughly the same uncertainty, as expected, and the weighted av- 
erage is 0.229 ± 0.004 Z , i.e. w 8% lower than the input value. 
This simple solution allowed us to correct for most of the under- 
estimate; however, a small (i.e. « 0.02 Z ) systematic still affects 
our measurements, especially in the outer regions. We also per- 
formed more realistic simulations, by introducing a background, 
and obtained substantially similar results. 

Although allowing observables to assume unphysical val- 
ues is against common sense, there are measurement procedures 
that can yield unphysical values. Deciding to accept only phys- 
ical values and reject others will clearly result in a bias. This 
is a general issue, which does not pertain to astrophysics alone. 
An interesting example we found concerns the analysis of data 
from the Collider Detector at Fermilat0. A statistical commit- 
tee specifically appointed to provide guidelines for the analysis 
of the Collider Detector data recommends a treatment similar to 
the one we propose here, i.e. unphysical values can be used in 
statistical procedures. 

The boundary effect we have just pointed out should play an 
important role when comparing subsamples characterized by a 
different statistical quality of the data. For example, in Sect. 14.21 
we compare the metallicity obtained for near and distant clus- 
ters. In our sample, near cluster spectra usually have a better sta- 
tistical quality for various reasons (e.g. longer observations, cos- 
mological dimming effect). If the metallicity is allowed to vary 
between and 5 Z Q , the mean metallicity beyond 0.20 R\so is 
0.249 + 0.011 Z Q and 0.188 + 0.014 Z , for near and distant clus- 
ters respectively. The measured discrepancy of 0.061 + 0.018 Z 
has a significance of more than three-sigma. Conversely, if the 
metallicity is allowed to vary between +5 Z (see Sect. 14.2b . the 
discrepancy of 0.029 ± 0.019 Z is consistent with a purely sta- 
tistical fluctuation. 

We then warn X-ray astronomers about the existence of this 
kind of systematics, which could affect the measurement of the 
metallicity, especially for the case of low metallicity, statistically 
poor spectra. 



with standard statistical quality). The curve is peaked on the in- 
put value, but the variance, which mainly depends on the number 
of counts around the energy of the emission line, is higher. Since 
the metallicity is forced to be positive, the curve is truncated 
and a dozen points pile up near the boundary (i.e. zero). In this 
case, the mean and the median are close to 0.25 Z , instead the 
weighted average is 0.193 + 0.003 Z , i.e. as 25% lower than the 
input value. Part of this discrepancy is due to a boundary effect. 
For measurements characterized by a negative fluctuation of the 
metallicity, the minimum of the^ 2 curve lies outside the allowed 
range. For these cases, the best fit values are forced to be above 
zero and the uncertainties have unreliable small values; there- 
fore, when computing a weighted average, these measurements 
have the highest statistical weight, and the net effect is an un- 
derestimate of the real metallicity. When fitting exactly the same 
spectra, allowing the metallicity to be negative, we obtain the 

histogram reported in the right panel of Fig. IA.2I For this case, 

the tails of the curve are more symmetric, all measurements have 5 http://www-cdf.fnal.gOv/physics/statistics/statistics_faq.html#ssell 



